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Abstract 

We propose a new mechanism for pattern formation based on the global al- 
ternation of two dynamics neither of which exhibits patterns. When driven by 
either one of the separate dynamics, the system goes to a spatially homoge- 
neous state associated with that dynamics. However, when the two dynamics 
are globally alternated sufficiently rapidly, the system exhibits stationary spa- 
tial patterns. Somewhat slower switching leads to oscillatory patterns. We 
support our findings by numerical simulations and discuss the results in terms 
of the symmetries of the system and the ratio of two relevant characteristic 
times, the switching period and the relaxation time to a homogeneous state 
in each separate dynamics. 
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Patterns in non-equilibrium systems arise in a variety of ways. For example, so-called 
dissipative structures arise from perturbations of homogeneous systems if the input of en- 
ergy is properly balanced by dissipation [jT| Perhaps the best-known dissipative structure is 
the Benard instability, formed when a layer of liquid is heated from below. At a given tem- 
perature, heat conduction starts to occur predominantly through convection, and regularly 
spaced, hexagonal convection cells are formed in the layer of liquid. This structure is present 
only as long as there is a supply of heat and disappears when this ceases. Many other exam- 
ples involve chemical oscillations in dissipative open systems [|l| . Patterns may also form by 
the temporal modulation of a parameter in systems that undergo Hopf bifurcations. This 
modulation may stabilitze standing waves or wave patterns that are otherwise unstable for 
any constant value of the parameter in this regime 0,^. A different sort of phenomenon, 
that of noise-induced pattern formation, may occur when an external spatially and tempo- 
rally uncorrelated noise affects a system parameter. The parameter fluctuations may induce 
spatial organization and pattern formation in a system where only a homogeneous state is 
stable when the parameter is constant [^^. 

In this Letter we identify and illustrate a new mechanism for pattern formation induced 
by a global periodic alternation between two dynamics, each of which by itself leads to a (dif- 
ferent) spatially homogeneous state. This mechanism is inspired by a number of examples 
where simple periodic switching between two dynamics, each of which produces "uninter- 
esting" or "disordered" or even "undesirable" states, leads to an "interesting" or "ordered" 
or "desirable" outcome. One example is the flashing Brownian ratchet, a rectifier of ther- 
mal fluctuations that can induce directed motion of Brownian particles merely by turning 
the ratchet potential on and off periodically 0. Another (in turn inspired by the ratchet 
example) is that of so-called paradoxical games, where the alternation of two losing games 
may lead to a winning game [^. Our proposed mechanism applies these ideas to spatially 
extended systems. 

We focus on a class of model systems based on the Swift-Hohenberg (SH) equation 
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^ (r, t) = ~V' (r, t)) + (r, t) + ^ (r, t) , (1) 

where the operator £ = — (1 + V^)^isa spatial couphng, and the white zero-centered Gaus- 
sian noise ^ (r, t) with correlation (r, t) (r', t')) = a^6{t — t')6{r — r') accounts for thermal 
or other fluctuations. The dynamics of the system is determined by the local potential V{ip), 
which in the standard SH model is an even quartic function of (p. The SH model leads to 
pattern formation (e.g. to the appearance of Rayleigh-Benard convective rolls) when the 
local potential has two stable equilibrium points. However, if V{(p) is monostable, no spatial 
structures appear in this system and the steady state is spatially homogeneous. The sta- 
bility boundary is identified by determining a uniform solution of the (noiseless) evolution 
equation (|l|) and linearizing about this solution. The stationarity condition V'{(f) + (f = 
(the term (p arises from the 1 in the coupling operator) has the solution Lp = (p. Setting 
(p = (p + Ap gives for the Fourier transform (denoted by a hat) of the linearized equation 

A^(k, t) = - [V"{^) + (1 - |k|2)2] A^(k, t), (2) 

and if V"{(p) < a morphological instability first sets in for the modes with |k| = 1. Note 
that in the usual SH model the even local potential leads to = 0, and that in general 
the patterns that emerge depend on the symmetry of the potential. Thus, symmetry of 
the evolution equation under the p < — > —p inversion leads to roll-shaped patterns, while 
absence of this symmetry leads to hexagonally-arranged spot-like structures 0. 

Consider now a global periodic switching of period T between two local potentials, Vi 
and V2. Thus, every point is subject to the same local potential at any given time. The 
system is governed by the equation 

V (r, t) = -vi (ip (r, t)) - fi (t) vi {p (r, t)) + Cp (r, t) + ^ (r, t) , (3) 

where V± {p) = (Vi (p) ± V2 (p)) /2 and fi (t) takes on the value 1 if t mod T < T/2 and — 1 
if t mod T > T/2. The resulting equation is a SH model with dichotomous periodic external 
forcing. To the best of our knowledge, this forcing differs from those previously considered 
in the literature 
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In order to investigate whether the switching mechanism can lead to pattern formation, 
we focus on local potentials Vi{ip) and V2{(p) that satisfy the stability conditions, 

V^{^i) + ^i = and V^'{^i)>0, (4) 

(monostable potentials) whereas their average V-^{(f) does not. Note that the condition 
that V+{(f) not satisfy (H) implies that both local potentials can not be quadratic. As a 
consequence, for the proposed pattern formation mechanism non-linearity is a necessary 
component. 

When each potential acts separately, the system tends to a homogeneous state. The 
expectation of pattern formation when the potentials are periodically alternated arises as 
follows. If the switching period T is large, every point in the system has time to equilibrate to 
the local potential before the local potentials switch. Regardless of the initial distribution, 
the entire system is expected to oscillate with period T between the homogeneous states 
corresponding to each potential. On the other hand, if T is sufficiently small (see below), 
^{t) can be adiabatically eliminated and replaced by its average value, which is zero. The 
system is then driven by the potential V^{(p), for which there is at least one state (f with 

Vl{ip) + ^ = and V"{ip) < 0, (5) 

and thus patterns appropriate to this local potential are expected to occur. Away from these 
extremes, the behavior depends on the switching rate. More specifically, the crossover time 
tr between slow and fast switching is the smaller of ti^2 and t2^i, where tj^j is the relaxation 
time, under the action ofVj, of the homogeneous state associated with V^,. We can estimate 
ti^j by focusing only on the k = mode and assuming that, when the potential switches 
from Vi to Vj, the mode amplitude behaves as a Brownian particle initially equilibrated in 
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the effective local potential Vi (</?) = Vi (</?) + When the local potential is switched, this 
point, which up to that moment was stable, becomes unstable. The relaxation time to the 
new homogeneous state associated with Vj is the time that it takes the Brownian particle 



to roll down the potential hill to the new equilibrium point |jTO| 



^*-i = ^/ dyexp d2;exp . (6) 

The behavior of the system (P) can be characterized by the ratio of the time T/2 that the 
system spends in each dynamics to the crossover time tr, r = T/2tr. 

Let us focus on a particular choice of potentials satisfying the conditions (|) and (j^): 

V.A^)^A,.{^±^-^T^), (7) 

where Ai 2 are positive constants. The corresponding effective local potentials with Ai = 
A2 = 1 are shown in Fig. |l[ The average potential is 

yj^ 09^ 09^ 

V+ {ip) = a+— + a-— - - (^-f (8) 

where a± = {Ai ± A2) /2. If r ^ 1 the entire system alternates between the homogeneous 
states (^1^2 • When r < 1, where we expect pattern formation with outcomes dependent on 
the specific value of r and on the equality or inequality of the parameters Ai^2- We support 
our reasoning with 2d numerical simulations of a discretized form of Eq. (|^) with these 
potentials, periodic boundary conditions, and parameters L^ = Ly = 64, a = 10~^ (the 
fluctuations must be sufficiently small not to swamp the potential barrier in V"+). Since the 
most unstable mode is |k*| ~ 1 [^, the typical wavelength of any pattern is expected to be 
A = 27t/ |k*| ~ 27r and the aspect ratio L/A ~ 10. In our simulations we either take the 
initial field to be random according to a Gaussian distribution (in which case the additive 
fiuctuations can actually be omitted entirely), or we can take an arbitrary initial condition 
(e.g. all points equilibrated with Vi), in which case the fiuctuations will distribute the field 
in any case. In some cases we choose an initial configuration that facilitates arrival at a 
particular final state simply to avoid a very long simulation time. At the initial time every 
local potential is set to, say, Vi, at time T/2 all the potentials are switched to V2, and so on. 
If Ai = A2, {(f) is an even function and, in the limit r ^ 0, Eq. (|^) satisfies inversion 

symmetry under the transformation (f < > —(f. In this case we expect the appearance of 

stationary rolls This is clearly seen in Fig. ^. Note that the widths of the rolls are 



consistent with the aspect ratio given earher. On an extremely long time scale, as in the SH 
model, the rolls line up in a more ordered fashion. 

As r increases toward the "resonance" condition r « 1, the contribution of {(f) can 
no longer be neglected. Hence Eq. (|^) will lack the symmetry (/? < — > —ip, and an oscil- 
latory spot-like pattern is expected 0. Furthermore, with Ai = A2, Eq. is invariant 
under the combined transformation {ip < > —'P,fJ' * ^ ~/^}) which requires a square spa- 
tial arrangement of the oscillatory pattern P|. A realization is shown in Fig. for which 
tr = ti^2 = ^2^1 ~ 2. The resonant period of the forcing is T ^ 4. We show color encoded 
snapshots of the field, where the oscillatory square patterns are clearly visible. Note that 
the field oscillates between a square lattice and its glide-transformed one. In other words, 
the spot centers do not move or oscillate; it is the surrounding background that oscillates. 
We are able to reproduce this behavior analytically through a decomposition in a small 
number of modes [|lT|. The size of the pattern units is again consistent with the aspect ratio 
given earlier. The oscillatory pattern and the lattice arrangement of the spots resemble the 
so-called oscillons found in vibrating granular media and clay [ |12| . 

The behavior is in some ways simpler when Ai ^ A2. The inversion symmetry ^ < > —(f 

is no longer satisfied. Therefore, both stationary (r <^ 1) and dynamic (r ^ 1) hexagonal 
spot patterns are expected [^. We have simulated the case Ai = 1, A2 = 2, for which 
the relaxation time is calculated to be ~ 1.6. As always, if r ^ 1 no spatial structures 
develop; homogeneous states simply alternate in time. On the other hand, for r = 0.94, we 
again obtain oscillatory patterns, as shown in the snapshots of Fig. ^ It is worth noting 
that this excitation density is quite different from the one in Fig. |^. As expected, the spots 
are arranged hexagonally. Most strikingly, in this case there is no glide oscillation but rather 
a true oscillation of localized excitations whose size is again consistent with the calculated 
aspect ratio. As r is decreased the excitation field is frozen and one obtains a stationary 
hexagonal pattern of spots. 

We have also performed simulations of Eq. when the modulation fi{t) is a dichotomous 
noise. In this case, the intermediate regime of oscillatory patterns fades out and the system 
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exhibits alternating homogeneous states if the correlation time of is large, and stable 
patterns (rolls or hexagons) if the correlation time is small |Tl|. 

Summarizing, we have shown that the alternation of two dynamics can create patterns 
even though each separate dynamics drives the system to a homogeneous state. Moreover, 
we have seen that in the crossover between alternating homogeneous states and stationary 
patterns, a rich phenomenology of oscillatory patterns may emerge. The mechanism is very 
general and can easily be extended to other situations such as reaction-diffusion systems. 
Thus, we have provided an alternative non-equilibrium mechanism to dissipative structures 
or modulation of bifurcation parameters for the formation of stationary and oscillatory 
spatial patterns. 

An interesting open problem is whether it is possible to find this type of behavior when 
the switching occurs between two equilibrium dynamics. This is not the case in the SH equa- 
tion, even when the homogeneous state is stable. For instance, in the case of convection, the 
homogeneous state describes a fluid in mechanical but not in thermal equilibrium. If switch- 
ing between local potential is achievable in the case of convection, our results would imply 
that the system will exhibit convection patterns by alternating dynamics which separately 
drive the fluid to mechanical (but not to thermal) equilibrium. 

The authors wish to thank R. Kawai for fruitful discussions. This work was partially sup- 
ported by the National Science Foundation under grant No. PHY-9970699, by DGES-Spain 
Grant PB-97-0076, by the New Del Amo Program, and by MECD-Spain Grant EX2001- 
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FIG. 1. Effective local potentials V\{}p) (solid curve) and V2(y') (dotted curve) witfi 
Ai = A2 = 1. Tfie mirror symmetry is broken if Ai ^ A2. 
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FIG. 2. Density plot of the field for the case Ai = A2 = 1 and r = 0.25. The roll-shaped 
pattern is essentially (see text) stationary. 
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FIG. 3. Snapshots of the field during a full period of the dichotomous forcing (^i = A2 = 1 and 
r = 1.15). Note the oscillations of the field between the square lattice and its glide-transformed 
one. 
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FIG. 4. Snapshots of the field for the case Ai = 1, A2 = 2 and r = 0.94 during a full period of 
the forcing function. The localized excitations are arranged in a hexagonal lattice. 
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